VORTEX LATTICE MELTING AND THE DAMPING OF THE DHVA 
OSCILLATIONS IN THE MIXED STATE 
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Phase fluctuations in the superconducting order param- 
eter, which are responsible for the melting of the Abrikosov 
vortex lattice below the mean field H C 2 , are shown to dra- 
matically enhance the scattering of quasi-particles by the fluc- 
tuating pair potential , thus leading to enhanced damping of 
the dHvA oscillations in the liquid mixed state. This effect is 
shown to quantitatively account for the detailed field depen- 
dence of the dHvA amplitude observed recently in the mixed 
state of a Quasi 2D organic SC. 

PACS numbers: 74.70.Kn, 74.60.-w, 74.40. +k 



Much effort has been recently invested in elucidating 
the fundamental mechanisms in which the Abrikosov vor- 
tex lattice is involved in the damping process of the de- 
Haas van-Alphen (dHvA) oscillations in the mixed state 
of pure, type-II superconductors. Several theoretical at- 
tempts have been made to account quantitatively for 
this effect (see discussions in Refs.EJJa ). It seems , how- 
ever, that non of the proposed theories is capable of fully 
achieving this goal, in particular, since there has been 
growing evidence recently to the significant sensitivity of 
the dHvA amplitude to extrinsic factors, such as vortex 
lattice disorder and flux lines sinningajd. Furthermore, 
it has been shown theoretically! 2 ] that any mechanism de- 
stroying the phase coherence in the vortex lattice, such as 
disorder or vortex lines fluctuations, should dramatically 
enhance the inhomogeneous Landau level broadening of 
quasi particles. 

The effect of vortex lattice melting is of special in- 
terest since it is associated with phase fluctuations of 
the superconducting order parameter, which can be im- 
plemented selfconsistently into .the Gorkov-Ginzburg- 
Landau (GGL) formalism used irJj,cl. 

In a type-II superconductor under a magnetic field a 
soft shear Goldstone mode, which .-can be described by 
long wavelength phase fiuctuationsa is connected to the 
Abrikosov lattice melting. Unfortunately rigorous ana- 
lytical approaches to this problem have encountered fun- 
damental difficulties: large order high temperature per- 
turbation expansion with Borel-Pade approximants to 
the low temperature behaviorEl,Q has no indication of an 
ordered vortex lattice even at zero temperature. The ex- 
isting non-perturbative approaches have not completely, 
clarified the situation: Rcnormalization group studiesu 



have predicted no crystal vortex state in a pure 2D su- 
perconductor (SC) at finite temperature, while tbp_novcl 
functional integral formalism suggested in Ref.EljLj, has 
led to some kind of a vortex liquid freezing transition 
without breaking the £7(1) symmetry. ■— Several Monte 
Carlo simulations, have recently showr£irli3 that even 
in a 2D SC a true vortex lattice melting phase transition 
takes place at finite temperature and that the transition 
is of the first order. 

In the present letter we make the first attempt to take 
into account the effect of vortex lattice melting on the 
magnetization (dHvA) oscillations in the mixed state. A 
simple, non-perturbative analytical model of melting is 
developed for this purpose within the GGL theory for a 
pure 2D extremely type II SC. Our model is based on 
the observation that in the vortex lattice state the main 
correction to the mean field free energy arises from fluc- 
tuating "Bragg chains", that is, from fluctuations which 
preserve long range periodic order along some crystallo- 
graphic direction in the vortex lattice. This simplification 
reduces our 2D fluctuation problem to a ID one, which 
is then solved exactly. 

Our starting point is the microscopic BCS Hamilto- 
nian for electrons interacting via an effective two-body 
attractive potential. We then write down the formal 
functional integral expression for the partition function of 
this system, eliminating the electronic field by introduc- 
ing bosonic Hubbard-Stratonovich complex field A (T*) 
, which describes all possible configurations of Cooper- 
pairs. By expanding the resulting log determinant in 
the small A(~f i '), up to fourth order, we then recover 
a rather general expression for the partition function in 
terms of the fully nonlocal Gorkov free energy functional 
F G [{A(r>),A*(rO}]@j!lnear H c2 : 



Z = J DA (T*) DA* (r > )exp [-F G /k B T] 



(1) 



The stationary phase approximation for this functional 
integral yields the mean field (GGL) equation for A 
In the lowest Landau level approximation this equa- 
tion can be solved by expanding in a set of y/~N Lan- 
dau functions (N being the total number of vortices) 
(f> n (x,y) = exp [iq n x - (y + q n /^) 2 } , with q n = a 
wave number along the x-axis, which is also proportional 
to the y-axis position of a given x -axis chain. Thus in 
the symmetric gauge 



1 



A(x,y)=e^ 



VN/2-1 

^2 a n<t>n (X, y) 
n=-VN/2 



(2) 



where for an arbitrary 2D periodic lattice the coefficients 

a „=<f!A exp (f 7 n*) I. All spatial variables are ex- 

pressed in units of the magnetic length. The parameter 
|A | 2 stands for the spatial average of the order param- 
eter square, and a x is the period along the £-axis. The 
summation in Eqj^ is performed over ^/N wavenumbers 
q n only. The reduction of the number of degrees of free- 
dom results from the assumed periodicity along the x 
-axis. Obviously a x is not unique and depends on the 
choice of the coordinate system; any 2D lattice can be 
described by two different values of a x , corresponding to 
the principal diagonals of the elementary cell (Fig.f). 

In Refsflp it was shown that the free energy functional 
near stationary solutions can be approximated by a lo- 
cal expression. This important property , resulting from 
gross cancellations of nonlocal terms in the free energy 
due to destructive interference between phase factors of 
the order parameter , simplifies the problem considerably 
, reducing it to the Ginsburg-Landau model with the well 
known local free energy functional 



F, 



GL 



-a|A(T>)| 2 + i/3|A(T>)| 4 



(3) 



Here a and (3 are known function of temperature and 
magnetic field (see below). It should be stressed, how- 
ever, that the locality of the free energy functional can 
be very sensitive to phase fluctuations, especially in the 
vortex liquid state , since these fluctuations tend to de- 
stroy the phase coherence of the order parameter. For 
the sake of simplicity we ignore in this paper possible 
nonlocality in the Boltzman factor of Eq.(Q). These as- 
sumption, though inconsistent with the nonlocal form of 
the free energy in the vortex liquid state, can be shown 
by a more careful analysis to have no significant influence 
on the final result. 

It is convenient to express the GL free energy through 
the Abrikosov parameter (3 a : Fgl — —Neo//3 a , where 
£o = net 2 / 2(3. In the lattice state (3 a is a function of a x 
and 7. The dependence of (3 a on z = 7r 2 /a 2 at 7 = tt/2 is 
plotted in Fig. (2a). Both minima with z\ = 7r/2v3 and 
Z2 = V3tt/2 correspond to the triangular lattice with 
(3 a ~ 1, 1596. The_maximum at 2:3 = tt/2 is obtained for 
the square latticeEil. 

In discussing the effect of fluctuations we first note that 
the order parameter in Eq.^ with arbitrary coefficients a n 
includes only those fluctuations which preserve periodic 
order along the x-direction. The number of independent 
coefficients a n is y/~N whereas the total number of orbital 
centers is N. Therefore each a n = \a n \ e i<Pn describes a 
set of \N orbital centers periodically arranged within 
a certain chain along the x-axis ("vortex Bragg chain"). 
The phase cp n determines the relative position x n = — 
Vn/Qn of the n-th chain. 



The scale of the melting temperature is determined 
by the value of the phase dependent terms in EqJ^. 
They can be readily calculated if we note that in 
the GL Hamiltonian there is a small parameter A = 
exp (— z), with z > 1 in the important regions near 
the minima of (3 a (see Fig. (2a)). Thus the quar- 
tic term can be written as an expansion in A, i.e. 



s'+f 



a*a„+ s a„ +t , while the leading phase 



dependent term is of the order A 2 . Integrating over am- 
plitude fluctuation \a n \ by using the stationary phase ap- 
proximation we find that the free energy can be writ- 
ten in the mean field like form (Fgl) — —Nso/(3 a , where 



A. 



1 



4A 2 (cos( Xn )) phas 



Xn 



4A 

tp n -i — 2ip n , and (■■■} phase means average over phase 

fluctuations. Note that fi = (cos (Xn)) phase = /o(tZ/t] ' 
where /i is the shear modulus of the Abrikosov lattice 
normalized by its mean field value. Here Ik (x) is a mod- 
ified Bessel function of the fc-order, and T cm = n^IXp 
is a characteristic crossover temperature from the crystal 
to the liquid state, with (3 m = (1 + 4A).Well below 



T 



1, so that thermal fluctuations do not distort 



significantly the triangular lattice and the Abrikosov pa- 
rameter is very close to its minimal (mean field) value, (3a, 
i.e. (3 a ~ (3 A ~ (1 + 4A - 4A 2 ). At T » T cm , on 
the other hand , [i — > and the vortex system trans- 
forms to a new (" liquid" ) state for which (3 a — (3 m - The 
crossover temperature T cm and the corresponding free 
energy depend on the choice of vortex chains (Fig.l). 
At zero temperature (i.e. in the vortex crystal state) 
the two minima at points z\ and zi have the same en- 
ergy (Fig. (2a)). At a finite temperature this degener- 
acy is removed by the phase dependent terms, which 
introduce different interchain coupling within different 
families of Bragg chains (corresponding to z\ and £2). 
Since z\ < z 2 , we have T cm (zi) > T cm (z 2 ). In Fig. (2b) 
we plot local minimal values of the average free energy 

(Fgl) as functions of the parameter t = — y^f?-, which 

has been used in RefR It is clear that at low tem- 
peratures the first state ( z\ ) is more stable than the 
second one ( Zi )■ Since T cm {z\) > T crn (22) , the free 
energy of the first state increases faster with increasing 
temperature than the second one, and so there is an in- 
tersection point T m at which energies of both states are 
equal but the corresponding entropies are a little differ- 
ent. Therefore we conclude that at this point there is 
a weak first order transition characterized by a small 
jump of the lattice entropy. The position of the cross- 
ing point t m ~ — 16 and the jump of entropy S = —T^ 
at T = T m : TAS/F MF ~ 7.5 x 1(T 3 , are in good agree- 
ment with the Monte-Carlo simulations^]. 

At higher temperatures (3 a does not describe cor- 
rectly the free energy. A straightforward calculation in 
this case yields for the partition function per unit vor- 
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/3 e 



erf (x) exp 



/ 3 In erf (a) y 



, where 



x = sign (a) v^j— . In contrast to mean field theory , in 
which there is a second order phase transition, in our case 
there is only a crossover to the normal metal state, with 
significant fluctuations of the superfluid density appear- 
ing at magnetic fields far above the mean field H C 2. 

All current theories of the dHvA effect in the mixed 
state are restricted to the mean field approximation. The 
situation in the vortex liquid state-is in some sense sim- 
ilar to that described by StephenEj in his derivation of 
Maki's damping formula. However, in this approach a 
random distribution of vortex lines was assumed, while 
the (nonzero) mean field value of the Abrikosov lattice 
order parameter was used in the calculation. The aver- 
aging over all realizations of the random vortex lattice in 
such a model yields for the quartic term of the free en- 
ergy: (a q+s+t a*a q+s a q+t ) ~ (S 8>0 + 5 t ,o). A similar ex- 
pression, which reflects the complete destruction of phase 
coherence in the vortex lattice, can be naturally and con- 
sistently derived in the model discussed here at T ^> T m 
by averaging over phase fluctuations, while at T > T m it 
is approximately valid. 

Thus, similarly to the MS theory, our fluctua- 
tion theory predicts an exponential damping: M sc = 



M n exp 



T 3/2 
Tj2 



A, 



where M sc , M n are ampli- 



tudes of the magnetization oscillations in the supercon- 
ducting and normal states respectively, Ao = Ao/huj c , 
and tif = Ef/?ilo c . A simple analytic expression for the 
damping parameter can be derived in the liquid state well 
above the melting point by integrating over amplitude 
fluctuations ( after neglecting small linear corrections in 
A): 



2 \ _ 1 gin Z v ^a/ d In erf (x) 
o| /7r da ~ 13 \ 2xdx 



(4) 



This expression has no singularity at the mean field 
transition, where both a — > and x — > , and smoothly 
interpolates between the high field (x —> —oo) value, 
T/tt I a I , and the low field (x — > oo) mean field 
value, a /(3. We have compared this damping with 
the experimental results on the organic quasi 2D SC: 
(ET) 2 Cu (SCN) 2 . In this computation the coefficients 
a and (3 derived previously in RefsD,0 for the GGL ex- 
pansion have been used, namely a = .5nhuj c In ^H^/H, 
and (3 = l.38nhuj c /nF , together with the well docu- 
mented Fermi surface and SC parameters of the studied 
material. The only adjustable parameter in this scheme 
is the mean field H C 2. We have found that a single value 
of about 4.7 T (see Fig. (3)) fits well the two different 
sets of available experimental, data, which were taken at 
quite different temperatureaijJiZl. 

It should be stressed that in the fluctuation theory the 
magnetic oscillations are smoothly damped well above 
H C 2 , and disappear below the mean field H C 2 in remark- 
able quantitative agreement with experiment. This con- 



trasts the result of the mean field MS theory , where the 
additional damping begins abruptly at H C 2 and then in- 
creases below H C 2 with significantly stronger rate than 
what is observed experimentally, 

The interesting point in our derivation of the four- 
particle correlation function concerns the fact that its 
factorization in the high temperature (or high field) limit 
into a product of two-particle correlation functions is 
solely due to integration over phase fluctuations, which 
leads to the limiting behavior:(a* +s+t a*a g+s a g+t )^ ?iase ~ 

(exp I (tp q+s + Lp q+t -<p q - Vq+s+t)) phase -> 

(S Si o + S t ,o) ■ 

This is a generalized form of the condition for the van- 
ishing of the shear modulus u in the liquid state. The 
apparent correlation between the vanishing of the shear 
modulus of the Abrikosov lattice and the factorization 
of the average phase factor of the quartic term in the 
SC free energy is a general result, independent of the 
nature of the melting transition. It is thus conceivable 
that at a first order melting transition, the vanishing of 
\x at some finite value of the parameter T/T m , m pre- 
dicted by the Monte-Carlo simulation reportedjnll3 ,will 
be accompanied by a transition from 'coherent'cl (i.e. rel- 
atively weak) damping in the crystal state to 'incoherent' 
MS-like damping in the liquid state at a certain magnetic 
field below H c2 . 

Such a transition will be difficult to observe in the quasi 
2D organic SC investigated due to the very small value 
of the melting temperature T m and to the strength of the 
2D fluctuations. In 3D superconductors where the role of 
fluctuations is less important and the melting transition 
is shifted significantly closer to H C 2 , one could expect 
such a transition to be observable. It is thus interesting 
to note that in the torque dHvA measurement performed 
on the borocarbide SC FA^-E^C" , where vortex lines 
pinning seemed to be weak (no peak effect near H C 2 ) , 
the Dingle plot exhibited a rather sharp upwarxLturn of 
its (negative) slope at a certain field below H C 2u3, which 
could indicate a freezing transition into an ordered vortex 
lattice stateEj. 

In conclusion, we have developed a simple analytical 
model of the vortex lattice melting in a 2D type-II SC, 
which is in good quantitative agreement with the state 
of the art numerical simulations in this field. We have 
shown that fluctuations in the SC order parameter, which 
are responsible for this melting, destroy the phase coher- 
ence in the quasi particle scattering and thus lead to en- 
hanced damping of the dHvA oscillations in the mixed 
state, in remarkable quantitative agreement with the ex- 
periment. 
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Figure Captions 
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Fig.(l): Two different choices of the vortex chains in 
the triangular lattice: dash lines - zi = ir/2^/3, solid lines 
- z 2 = y/Zit/2. 

Fig. (2): (a) Mean field Abrikosov parameter, /3 Q , as a 
function of z — (ir/a x ) 2 . (b) Dependence of the local 
minimal values of the free energy — (F) / Fmf >on the pa- 
rameter t = — 2y / £o/T. 

Fig. (3): The magnetization damping factor R s = 
M sc /M n , as a function of magnetic field H in the quasi 
2D SC:(ET) 2 Cu(SCN) 2 . Solid lines are our theoreti- 
cal curves for T = 20m_R',and for the dHvA frequency 
F = 690T (1), and T = 120mK, F = 600T (2);stars and 
circles are—tha corresponding experimental data taken 
from Refs.EjjHl, respectively. The other parameters used 
in the calculations are:T c = 10AK (from which Ao 
is calculated using BCS weak coupling formula), and 
to* = 3.5m e . 
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